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Abstract 

The gravitational evolution of the genus of the density field in large-scale structure is 
analytically studied in a weakly nonlinear regime using second-order perturbation theory. 
Weakly nonlinear evolution produces asymmetry in the symmetric genus curve for Gaus- 
sian initial density field. The effect of smoothing the density field in perturbation theory 
on the genus curve is also evaluated and gives the dependence of the asymmetry of the 
genus curve on spectra of initial fluctuations. 
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1. INTRODUCTION 



Recently redshift surveys have been unveiling new detailed information of structures 
in the universe. To characterize the pattern of the structures, the topological analysis 
of the galaxy distribution provides a useful and intuitively clear way. Gott, Melott & 
Dickinson (1986) proposed to use the Euler characteristic of surfaces of constant density 
as a quantitative measure for the topology of large-scale structure. The genus, which 
is widely used in subsequent topological analyses, is defined by —1/2 times the Euler 
characteristic per unit volume. More intuitively, the genus corresponds to "the number of 
holes" of the surfaces minus "the number of isolated regions" surrounded by the surfaces 
per unit volume. The genus is a function of smoothing scales and the density threshold. 
The genus as a function of density threshold for a fixed smoothing scale is called the genus 
curve and is analyzed both in numerical simulations and in redshift surveys of galaxies by 
many authors (Gott, Weinberg & Melott 1987; Weinberg, Gott & Melott 1987; Melott, 
Weinberg & Gott 1988; Gott et al. 1989; Park & Gott 1991; Park, Gott & da Costa 1992; 
Weinberg & Cole 1992; Moore et al. 1992; Vogeley, Park, Geller, Huchra & Gott 1994; 
Rhoads, Gott & Postman 1994). 

The only analytical expression for the genus curve known so far is for Gaussian random 
density field (Adler 1981; Doroshkevich 1970; Bardeen, Bond, Kaiser & Szalay 1986; and 
Hamilton, Gott & Weinberg 1986), and is given by 

G ^-h^T-T e "" 2{1 -^ (L1) 

where v is the difference between density threshold and mean density in units of standard 
deviation of density and 

{k } = JP(k)d^k ' (L2) 

with P(k) being the power spectrum of the density fluctuation. Previous analyses mainly 
compared the observational genus with the random Gaussian prediction ( |1 . 1| ) . With suf- 
ficiently large smoothing scales, this comparison could tell us if initial density fluctuation 
is random Gaussian. With the finite smoothing scale of cosmological interest, the ef- 
fect of nonlinear gravitational evolution on the genus curves would be substantial. This 
nonlinear effect has been explored only by using iV-body numerical simulations so far. 
The main purpose of this Letter is to approach analytically this problem in the weakly 
nonlinear regime using a second-order perturbation theory. In the following, the general 
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formula of the genus curves for the field with weak non-Gaussianity is presented. This 
formula is incorporated in the second order perturbation theory. The smoothing effect on 
perturbation theory is also considered. 



2. THE GENUS CURVE FOR QUASI-GAUSSIAN RANDOM FIELD 

Random Gaussian fields are characterized by the fact that connected correlation func- 
tions, except the second-order correlation function £, all vanish. This is why the genus 
curve (eq. ||l.l|| ) for random Gaussian field is completely determined by P(k) which is a 



Fourier transform of £. In this section, correction terms for equation (|1.1|) in the pres- 
ence of weak non-Gaussianity are presented. The term "weak non-Gaussianity" is defined 
below. 

To simplify the notation, the following seven quantities for a non-Gaussian random 
field 6(x, y, z) with zero mean are denoted by (// = 1, . . . , 7): 

8_ l_cM ld5_ im ld*5_ 1(P5_ 1 d 2 5 
a ' cr dx ' a dy ' a dz ' a dx 2 ' a dy 2 ' a dxdy ' 



where a = y (S 2 ) is an rms of the field and the field is defined in Cartesian coordinates 
x, y, z. The field 5 is identified with density contrast p/p — 1 in astrophysical applications. 
The Euler characteristic per unit volume n x (y) of constant surfaces 5 = vo is given by 
(Adler 1981; Bardeen et al. 1986) 

n x {y) = (p^A,)) , (2.2) 

where 

p x (A^ = 5 D (At - u)6 D (A 2 )6 D (A 3 )\A 4 \(A 5 A 6 - A 2 ), (2.3) 

and 8b is a Dirac's delta-function. This expression is valid for general non-Gaussian 
random fields. Let us proceed to evaluating the expectation value (eq. ||2.2j| ) for weakly 
non-Gaussian field. As usual, we define a partition function Z{Jp) as a Fourier transform 
of a distribution function P(A fl ) of quantities A^. 



Z(J^) = r d 7 AP(A^)e W (i^JuAu) . 

J-co \ u ) 



(2.4) 



The cumulant expansion theorem (e.g., Ma 1985) states that In Z{J^) is a generating func- 
tion of connected correlation function W^}.^ N = (A^ ■ ■ ■ A^ N ) C (see Bertschinger 1992). 
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Then one obtains the inverse Fourier transform of equation ( |2.4| ) in the following useful 
form: 

W = «p E 1 ^ £ Vff?.^ gr W, (2.5) 

where 

PcKJ = 1 exp ("£4. (V^" 1 )^) > ( 2 - 6 ) 

J(27r)Met V M ' 

is a multivariate Gaussian distribution function characterized by a correlation matrix . 
In a weakly non-Gaussian case, the exponential function in equation ( |2.5| ) is expanded 
in Taylor series and equation (|2.2|) is expanded by higher-order correlations. In the fol- 
lowing;, we assume that ^ {A ° ~ 0{a N - 2 ). This relation is a very definition of "weak 
non-Gaussianity" in this letter and is a result of perturbation theory (Fry 1984; Goroff et 
al. 1986; Bernardeau 1992). Thus to the first order in a, equation (|2.2j ) reduces to 

"xM - MA)) a + I £ e ( S Q^ r ) G + 0(a% (2.7) 

where (• • -)q denotes averaging by multivariate Gaussian distribution (eq. |[2.6|| ). All the 
terms in r.h.s. of equation ( |2.7| ) can be evaluated by straightforward but tedious Gaussian 
integrals. Spatial homogeneity and isotropy simplify the final result as 

„2\ \ 3/2 



H,(i>) + o ( -tf 5 (z/) + —H 3 (v) + 3UH 1 (u) ) + 0(a 2 ' 



(2.8) 

where H n (v) = (-)V / 2 (d/dv) n e~ u I 2 are Hermite polynomials, and we have defined 
three quantities, 



a 



^ = -7T77^-7<^A5>, (2-9) 



2{k 2 )a 

The quantity S is usually called "skewness". The first term in square brackets of equa- 
tion Q2.8| ) corresponds to Gaussian contribution and the other terms correspond to non- 
Gaussian contribution to Euler number density. 

As an illustrative application of this result, we consider the case that correlation 
functions are given by hierarchical model. In hierarchical model, connected correlation 
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function of N-th order is modeled as a sum of N — 1 products of £, thus our previous 
assumption ~ 0(a N ~ 2 ) is satisfied. Specifically, third order correlation function 
C(x, y, z) = (5(x)5(y)5(z)) is related to £(|ac - y\) = (5(x)5(y)) by 

C(x 1 y 1 z) = Q\i{\x-ymy-z\)+i{\y-zmz-x\) + i{\z-x\)i{\x-y\)i (2.10) 



where Q is a constant (undetermined in this model). If this equation ( 2.10| ) is exact for 
some large smoothing scale such that £(0) = d 2 < 1, the quantities S, T and U reduces to 
3Q, 2Q and Q, respectively, and the genus curve with correction terms of non-Gaussianity 
is 



G (hi.), 



1 



V) 



4tt 



I - z/ 2 - ? (z/ -4z/ 3 4 3/ ') a 



+ 0(a 2 ), (2.11) 



Figure 1(a) plots the results for Qa = 0, 0.2, 0.4, 0.6. 

3. GRAVITATIONAL EVOLUTION OF THE GENUS CURVE IN SECOND 

ORDER PERTURBATION THEORY 



Gravitational nonlinear evolution give rise to S, T, U even from the initial Gaussian 
random density fluctuation which has vanishing S, T, U. We use second order pertur- 
bation theory of the nonrelativistic collisionless self-gravitating system in the fluid limit 
(e.g., Peebles 1980, §18) to compute S, T, U to lowest order in a. Here we present the 
results explicitly in the Einstein-de Sitter case, Q = 1, A = 0. The effects of parameters 
fl, A on our results are expected to be weak as discussed in the next section. Considering 
growing mode only, third order correlation function in Fourier space is given by (Fry 1984; 
Goroff et al. 1986) 



5(fe 1 )5(fe 2 )5(fe 2 



12 + (h + M fc i fc 2 + 4 / fci-fc 2 ' 

7 \k 2 kij k\k 2 7 \ k x k 2 

x(27r) 3 ^(fc 1 + fc 2 + fc 3 ), 



P(h)P(k 2 ) + eye. 

(3.1) 



where P(k) is a power spectrum of linear theory. Evaluating S, T, U (eq. |[2.9|| ) in Fourier 
space with equation ( |3.1| ) results in, 



34 

y 



82 
21' 



U 



54 
35 



(3.2) 



- 5 - 



This value of the skewness S was already given by Peebles (1980). Note that S, T, U are 
independent on the shape of the initial power spectrum. The genus curve predicted by 
perturbation theory is, therefore, 



> 2 > 



3/2 



> /' 17 5 47 3 4 

V - — v V V I cr 

21 21 5 



+ (3.3) 



which is plotted in Figure 1(b) for a = 0, 0.2, 0.4, 0.6. 

In fact, the observable curve is obtained by smoothing of density fluctuation while 
the above result is not for smoothed field. Let us evaluate the smoothing effect in the 
case that the smoothed density fluctuation with sufficiently large smoothing scale is well 
described by second order perturbation theory. Recently, Juszkiewicz, Bouchet & Colombi 
(1993) gave the smoothing effect on the skewness S for Gaussian and top-hat filter. In the 
following, we use the Gaussian filter 5 R (x) = (2itR 2 )~ ?J l 2 J d 3 y5(y) exp(— \x — y\ 2 /2R 2 ) 
which is usually adopted in smoothing noisy observational data. The quantities S, T, U 
for Gaussian smoothed field 5r are obtained similarly as in the unsmoothed case, and the 
result is 

3 

S = 28^(R) 151220 + 7/131 + 2/222] ' 

T = u J a , {R) [ 10/ 240 + 12/331 + 71x51 + 11/242 + 2/ 333 ] , (3.4) 

flk 2 )\~ 2 1 

U = \T ) 168^) [5/440 + 7/351 ~ 3/442 ~ 7/353 ~ 2/444] ' 
where we defined 

/* 00 r 00 c 1 

Imnr = dx dy d^ 2 ^ 2 ^ x m y n ^ P (x) P (y) . (3.5) 
Jo Jo J-l 



Wheb R — > 0, the dependence of equation (|3.4|) on the initial power spectrum is caceled 
and equation ( |3.2|) is rederived. Thus the genus curve for smoothed field is dependent 
on the specific shape of initial power spectrum on the contrary to the unsmoothed one. 
For the case of initial power spectrum with power-law, P(k) oc k n , equations ( |3.4| ) are 
expressed using generalized hypergeometric functions p F q as 

_ 30 n /n + 3 n + 3 3 1\ , o . ^ fn + 3 n + 5 5 I 

S = ^2^1 — , ^ t - ( n + 3 )2^i 



7 V 2 ' 2 7 2 7 A) y ' V 2 1 2 1 T 4 

4 /n + 3n + 3315 1 
+ -3/2 



T~^\ 2 ' 2 '2'2'2'4 
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m 20 „ /n + 3 n + 5 3 1\ 4 „ /n + 5 n + 5 5 1 
T = — 2 F l — — , — — ; -; - - ~ 2 iq 



7 x V 2 ' 2 ' 2' 4/ 7 V 2 ' 2 ' 2' 4, 
n + 5 / n + 3 n + 7 5 1\ 22 / n + 3 n + 5 3 15 1, 

3 V 2 ' 2 '2'4y + 21 3 2 V 2 ' 2 ' 2' 2' 2' 4 ' 
2(n + 3) /n + 5 n + 5 5 3 7 1 



35 " 3 2 V 2 ' 2 ' 2' 2' 2' 4 

15 /n + 5 n + 5 3 1\ n + 5 /n + 5 n + 7 5 1 
U = -772F1 — - — , — - — ; -; — 2F1 



7 V 2 2 ' 2' 4/ 2 V 2 2 24 

3 /71 + 5 n + 5 3 1 5 1\ 3(n + 5) /n + 5 n + 7 5 3 7 1 

~7 3 2 ' ~ 2~ ' 2' 2' 2' 4 J + 10 3 2 ' ~ 2~ ' 2' 2' 2' 4 
6 /n + 5 n + 5 5 1 7 1 

— — 3-^2 



35 V 2 2 '2 2 2 4, 
Table 1 shows the numerical values of these quantities for specific values of n. The 
corresponding genus curves for n = 1,0,-1,-2 are shown in Figure 1(c) to 1(f) for 
a = 0,0.2,0.4,0.6. 



4. DISCUSSION 

As seen in equation ( |2.8| ), the nonlinear correction of first order in a is an odd function 
of v and generates asymmetry between high-density region and low-density region in the 
genus curve, though it does not change the amplitude G(0). The pattern of the asymmetry 
in the genus curve of smoothed density field is dependent on initial power spectra. Thus, in 
principle, observations of the genus curve can restrict the properties of initial fluctuation, 
such as Gaussianity, the shape of the spectrum, by the amplitude and the pattern of 
asymmetry of the curve. The presently available redshift data of galaxies are not enough 
to have the statistically sufficient accuracy on topology of the large-scale structure. The 
projects as Digital Sky Survey (DSS), however, will enable us to have a large amount of 
redshift data in near future and the analysis indicated in this letter will be important one. 

We have shown the results explicitly in Einstein-de Sitter case. Bernardeau (1993) 
obtained skewness S smoothed by the top-hat filter for arbitrary Q and A and found that 
the dependence of skewness on these parameters is extremely weak. Applying his method, 
the Gaussian filtered S, T and U can also be evaluated for arbitrary Q and A. Explicit 
calculations show that the dependence on these cosmological parameters is also weak and 
the Einstein-de Sitter case is a good approximation to the other cases (we will give the 
explicit evaluations elsewhere). 

The effect of biasing (Kaiser 1984; Bardeen et al. 1986) between the galaxy distribution 
and the matter distribution on the genus curve is an important issue. Quite generally, 
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the correlation functions arose from local bias approach to the hierarchical model in the 
large-scale limit (Szalay 1988; Fry & Gaztanaga 1993; Matsubara 1994). The parameter 
Q in hierarchical model is determined by biasing mechanism and can take the positive and 
negative values. Thus the effect of biasing on the genus curve is approximately expressed 
by equation ( 2.11j ). 

The method used in this letter gives the general way to express the local function 
of density field by correlation functions: we can choose any function instead of equation 
(|2.3| ). For example, the evaluation of nonlinear effect on level crossing statistics (Ryden 
1988; Ryden et al. 1988) are straightforward using our method. In the next work, gen- 
eralizations to the nonlocal functions of density field and the multi-point statistics are 
elegantly described by diagrammatic language which provides a powerful way to the sta- 
tistical analyses in large-scale structure in the universe (Matsubara 1994). Applications 
of techniques developed in this letter to the statistics of anisotropy of cosmic microwave 
background radiation map are now under progress in search for non-Gaussianity of pri- 
mordial fluctuations in the universe. 



I am grateful to thank Y. Suto for a careful reading of the manuscript and useful 
comments. 
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unsmoothed 


n = 1 


n = 


n= -1 


n = -2 


n = — 3 


s 


4.857 


3.029 


3.144 


3.468 


4.022 


4.857 


T 


3.905 


2.020 


2.096 


2.312 


2.681 


3.238 


U 


1.543 


1.431 


1.292 


1.227 


1.222 


1.272 



Table 1: The numerical values of S, T, U for unsmoothed perturbation theory and 
smoothed perturbation theory for power-law spectra, n — 1 to —3. 
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FIGURE CAPTIONS 



Figure 1 : Asymmetries of the normalized genus curves G(v)/G(0) induced by non- 
Gaussianity to first order of rms a of fluctuation. The sources of non-Gaussianity 
are (a) hierarchical model, (b) unsmoothed perturbation theory with Gaussian initial 
fluctuations, Gaussian-smoothed perturbation theory with Gaussian initial fluctua- 
tions of power-law spectra, (c) P(k) oc k, (d) P(k) = const., (e) P(k) oc fc -1 , (f) 
P(k) oc k~ 2 . Solid lines, dotted lines, dashed lines, long-dashed lines show a = (in 
(a), Qa =) 0,0.2,0.4,0.6 cases, respectively. 
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Fig. 1 
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_ a) Hierarchical model 
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_ b) P.T. without smoothing 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
_ c) n = 1 with smoothing 











d) n = with smoothing 
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e) n = -1 with smoothing 
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f) n = -2 with smoothing 
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